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Testing the GRB Variability/Peak Luminosity correlation using the 
pseudo-redshifts of a large sample of BATSE GRBs 



in 
o 
o 

< 

m 
(N 



> 
cn 
oo 

00 

o 
in 
o 

o 



Cristiano Guidorzi^* 

^Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Birkenhead, CH41 ILD, UK 



1 INTRODUCTION 



Years after t he discovery of the firs t X-ray jCosta etalJll997h 
and optical I Van Paradiis et all Il997h afterglow counterparts to 
Gam ma-Ray Bursts (GRB s) and the determination of the first red- 
shift iMetzger et alj[T997l) . it has been possible to discover some 
correlations between burst-frame properties derived on a sample of 
some tens of GRBs with measured spectroscopic redshift. Some 
of them involve spectral and temporal properties: e.g., the anti- 
co rrelation between p eak luminosity and spectral lag discovered 
bv lNorris et alJ |200y), and the correla tion between peak luminos- 
ity and temporal variability found bv iFenimore & Ramirez-Rui3 
(2000) and Reichait et al. {200^ (hereafter FRROO and ROl, re- 
spectively). Other correlati ons concern the sp ectrum and the ener- 
getics of the GRBs, like the lAmati et alji2002l) correlation between 
the peak energy Ep of the burst-frame E F{E) s pectrum and the 
total isotropic-equivalent released energy, or the 'Ghirlanda e t all 
(2004) one between Ep and the collimation-corrected released en- 
ergy for those GRBs for which it was possible t o estimate the 
beaming angle. Similarly to the Amati relationship, Yonetoku et alj 
i2004) found that Ep also correlates with the peak luminosity L. 
Some of these relationships have been used to estimate the red- 
shifts of large samples of GRBs. In particular, .Band et al. (2004h 
(hereafter BNB04) used the anticorrelation between peak luminos- 
ity and spectral lag to estimate the redshift of 1 186 BATSE GRBs. 
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ABSTRACT 

We test the correlation found bv'Reic hart et alJ ('2001') between time variability and peak lu- 
minosity of Gamma-Ray Bursts (GRBs). Recently Guidorzi et al. (2005) found that this still 
holds for a sample of 3 2 GRBs with spectros copic reds hift, although with a la rger scatter than 
that originally found bv lReichart et alJ ( 1200 10 . However lGuidorzi et alJ (l2005l) also found that 
a power law does not provide a good description of that. We report on the same test per- 
formed on a sample of 551 BATSE GRBs with a significant measure of variability assuming 
the pseudo-redshifts derived by Band et al. (2004) (1186 GRBs) through the anticorrelation 
between s pectral lag and peak l uminosity. We still find a correlation between variability as 
defined by Reichart et al. ('2001') and peak luminosity with higher significance. However, this 
subsample of BATSE GRBs show a higher scatter around the best-fitting power law than that 
found by Reichart et al. (2001) in the variability/peak luminosity space. This is in agreement 
with the result found by Gi iidorzi et a l. ( 2005) on a sample of 32 GRBs with measured red- 
shift. These results confirm that a power law does not provide a sa t isfacto ry description for all 
the GRBs, in contrast with the original findings bv lReichart et alJ ( 1200 ll) . 
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Recently, the deba te on the Amati and G hirla nda correlations 
saw a couple of papers. iNakar & PiranI i20o'^ and IB and & Preec3 
i2005h . according to which at least a considerable fraction of the 
overall catalogue of BATSE GRBs cannot be consistent with those 
correlations. On the other side, three other papers suggest that 
the majority of the BATSE GRBs ar e cons istent with the Amati 
and G hirla nda correlations : Ghirland a et alJ 1200 5). Bosni ak et alj 
i2005h and IPizzichini et alji (,2005.) . These authors, starting from 
different assumptions concerning the redshifts of the BATSE 
GRBs, independently show that the Amati and Ghirlanda correla- 
tions are confirmed. See the review by Amati (2005) for a detailed 
discussion of the debate. 

Concerning the correlation originally found by ROl between 
variability and peak luminosity , als o confirmed fo r X-Ray Flashes 
(XRFs; see 'Heise et al.' bOOll)) bv lReichar t et al." (2003). recently 
lGuidorzi"et al. (2005) (hereafter G05) have confirmed it through a 
larger sample of 32 GRBs with measured redshift, although with 
several differences from ROl. In fact, they find a much larger scat- 
ter than that found by ROl with some notable outlier. Furthermore, 
G05 found that the p ower-law description originally obtained by 
iReichart et alj J200ll) for a sample of 13 GRBs is no more satisfac- 
tory for an enlarged sample of 32 GRBs with measured redshift. 

In this paper we test the variability vs. peak luminosity corre- 
lation similarly to G05, this time using a large sample of BATSE 
GRBs assuming the pseudo-redshifts deri ved from BNB04 . In par - 
ticular, we used a different approach from lGhirlanda et'al] i2005h : 
since BNB04 do not provide confidence intervals on both redshift z 
and peak luminosity L, following BNB04 we assumed the validity 
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of the lNorris etall J200dl anticorrelation between peak luminosity 
and burst-frame spectral lag and for each BNB04 GRB we derived 
confidence intervals on both z and L and made sure that they were 
consistent with the BNB04 catalogue. Eventually we compared the 
results with those by G05. 

The idea of using a large sample of BATSE GRBs with un- 
known reds hift to test the variability/peak luminosity correlation 
is not new: "Schaefer et al.' (2001) used 112 BATSE GRBs to test 
the correlation between spectral lag and variability, as defined by 
FRROO, as both have been shown to be correlated with peak lu- 
minosity. The motivation of the investigation reported in this pa- 
per is based on several aspects. Firstly, it must be pointed out 
that Scha efer et alJ ^00?) made use of different definitions of both 
variability and peak luminosity, very similar to those given by 
FRROO. Although the two definitions of variability given by ROl 
and FRROO appear to be correlated (see fig. 3 in FRROO), the rela- 
tion between the two does not seem to be direct proportionality. In 
fact. lSchaefer et alJi2001D found a best-fitting power-law index be- 
tween L and 1^ of msoi = 2.5± 1.0 to be compared with that found 
by ROl, mRoi = S.SlIiJ'g. Secondly, nowadays the greater num- 
ber of GRBs with known redshift than that available at that time, 
allowed to refine the best-fitting parameters of the lag/luminosity 
correlation tNorris 2002). For the first time, here we compare the 
properties in the variability/peak luminosity space of a sample of 32 
GRBs with measured redshift studied by G05 with those of a large 
BATSE GRBs sample, using the definitions by ROl for variability 
and peak luminosity. 

In Section|2|we discuss how we derived the samples of GRBs. 
In Section|3|we describe how we calculated the variability. In Sec- 
tion|4|we present our results and in Section|5|we discuss them. 



2 THE GRB SAMPLE 

From the catalogue published by BNB04 we first selected 1186 
GRBs, for which all the information derived by BNB04 is avail- 
able: redshift, spectral lag, best-fitting spectral parameters, peak 
luminosity. After requiring a significantly positive time lag, this 
sample shrank to 866 GRBs. Since BNB04 do not provide un- 
certainties on their measures of redshift and peak luminosity, we 
evaluated them as follows. We used the photon peak count rates 
(ph cm~'^ s~^) in the 50-300 keV band measured on a 256-ms 
time-scale, as reported in the BATSE catalogue ^ . In order to eval- 
uate the energy peak flux (erg cm~^ s~^) we need to know the 
spectrum. We adopted the same values as BNB04: i.e., we made 
use of the best-fitting spectral parameters found bv lMallozzi et alJ 
il99 8h who fitted the energy spectrum of the peak with the Band 
function (Band et al. 1993). For those GRBs for which this piece 
of information was not avail able, l i ke B NB04 we also assumed 
the average values found by Preecd i2000 ) a = —0.8 ± 0.1 and 
P = -2.3 ± 0.1, while E-p = Eq {2 + a) was taken from 
BNB04. From the photon peak flux we calculated the normalisa- 
tion. We then evaluated the bolometric energy peak flux Fb in the 
0-10* keV band and its uncertainty. 

We calculated the best value and a 2-a confidence interval for 
the redshift of each GRB assuming the validity of the anticorrela- 
tion between burst-frame spectral lag tb and bolometric peak lu- 
mi nosity Lb orig inally found bv iNorriret al. (.200 0) and refined 
bv lNorriji2003) : L50 = 21.8 (tb/O.SS s)"^ ". Lb = L50 x 



^ http://gammaray.msfc.nasa.gov/batse/grb/catalog/current/tables/flux_table.txt 
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Figure 1. Example of how we determined the redshift confidence interval 
in the case of trigger # 228. Dashed line shows the Noms relation. EiTor 



bars are at 2-(7. The result for this GRB is z 
for this GRB ^ = 1.435. 



1.77^ 
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lO^'' erg s~^. We determined the range for the redshift z for which 
the point {tb,Lb) was consistent within 2-a with the Norris rela- 
tion (see Fig.0. We used the following: re = ro/(l + 2:)'^ (ro is the 
spectral lag measured in the observer frame; c = 0.6, that accounts 
for both cosmol ofiical dilation and narr owing of pulses with higher 
energy found bv lFenimore et alj il995t) ) and Lb = 47r Df^ (2) _Fb, 
where Dl{z) is the luminosity distance at redshift z. Here we 
adopted the same cosmology as BNB04: Ho — 70 km s^^ Mpc^^, 
Qm = 0.3, and = 0.7. 

We point out that the uncertainties on tb and Lb have been 
propagated from those affecting the measures of the lag itself and of 
the energy peak flux, respectively (at each fixed z). We studied the 
scatter distribution of the 2-a lower and upper limits with respect 
to the Norris power-law: it turned out that for the 67% of the entire 
sample these limits are s cattered more than 0. 26 around the best- 
fitting power law. Since ISchaefer et alj i200lh found a scatter of 
0.26 we conclude that at least a fraction of 67% have estimated 
luminosities that are consistent with the scatter intrinsic to the lag- 
luminosity relation. 

This is differe nt from the approach followed by 
iGhirlanda et alj i2005l ). who assigned each GRB peak lumi- 
nosity an uncertainty derived from its scatter with respect to the 
Norris relation (see their eq. 5). Understandably, our confidence 
intervals for z are fully consistent with the values provided by 
BNB04 fo r the entire samp l e of 8 66 GRBs so far considered. 

Like I Ghirlanda eiei\ i2005. ). we excluded 38 GRBs with 
Ep ^ 40 keV because close to the lower boundary of BATSE en- 
ergy pass-band and therefore likely to be biased. In addition, one 
GRB was rejected for its undetermined redshift range, leaving us 
with a BNB04 final subset of 827 GRBs. 

Eventually, in order to comply with ROl energy band and cos- 
mology, making use of the same spectral parameters as above, we 
derived the best peak luminosity and its 2-a confidence interval in 
the burst-frame 100-1000 keV energy band from the best value for 
z and its 2-a confidence interval, assuming this time Hq — 65 km 
s~^ Mpc"^ like ROl and G05. Let be the photon spectrum at 
peak (ph cm^^ s^^); L in the burst-frame 100-1000 keV was then 
computed according to eq.Q 



L = A-KDliz) 



1000/(l-l-z) 



dE 



(1) 



100/(1 + 2) 
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This is formally the same as eq. 9 in ROl and eq. 8 in G05. 



3 VARIABILITY ESTIMATION 

Variability V was estimated using the public BATSE 64-ms con- 
catenated light curves ^. For each GRB we interpolated the back- 
ground by fitting with polynomials of up to 4"^ degree as prescribed 
by the BATSE team . This procedure was applied independently to 
each of the four energy channels: 25-55 keV, 55-1 10 keV, 1 10-320 
keV, and > 320 keV. Variability was computed for each GRB inde- 
pendently in each channel according to the definition by ROl (eqs. 
4-8 therein). Following ROl, we adopted a smoothing time-scale 
of Tf with / = 0.45: T/ is the shortest cumulative time interval in 
which a fraction / of the total fluence is collected. Hereafter Vf is 
the variability obtained adopting a time-scale of Tf. Uncertainties 
on Vf have been calculated combining the statistical uncertainty 
expressed by eq. 8 of ROl with that due to the error on redshift z. 

Eventually, for each GRB we performed a test and re- 
jected all the GRBs showing significantly different variability mea- 
sures between different energy channels. This requirement relies 
on the definition of Vf: as explained by ROl, the definition of 
variability already accounts for the narrowing of pulses with en- 
ergy I Fenimore et al. 1995) (pulses' width is proportional \.o E^" , 
a = 0.4) and for the cosmological energy shift. For those GRBs 
showing consistent measures of V/=o.45 in all channels we consid- 
ered the weighted average. 




Scatter ALog(L5o) 

Figure 4. Distribution of the scatter of the two groups around the best-fitting 
power law obtained by G05: 551 BNB04 GRBs (black solid line) and 31 
GRBs with measured redshift (red dashed line; the peculiar subluminous 
GRB 980425 is not considered here). Also shown are the two best-fitting 
normal distributions, whose cr's are 0.70 (BNB04) and 0.57 (G05). 



Table 1. Best-fitting parameters for the two GRB sets for which we find 
significant correlation between variability and peak luminosity, m and q 
are the parameters of the best-fitting power-law: L50 = IC VJ^q 45- The 
last column reports the scatter around the best-fitting power law. 



GRB Set m q -x^lAof scatter 



4 RESULTS 

The selection of the GRBs with a significant and consistent mea- 
sure of variability reduced the sample from 827 to 551 BNB04 
GRBs. 

Figure |2| plots V>=o.45 vs. L for the subset of 551 BNB04 
GRBs. We also plot 31 GRBs with known redshift derived by 
G05 using d a ta fro m GKBMI BeppoSAX (Ferocietal. 1997; 
iFrontera et all Il997^ BPSSEICGRO (Paciesas et al. 199^ 
FREG ATE/Z/gra:-// iAtteia et alj l20()3h . Ulysses ( .Hurlev et "Si 
Il992h. Konus/My P jlptekar et alJ Il995) and BJ^ I Swift 
jGehrels et all2004^ . From the sample of 32 GRBs studied by G05 
we ignored the peculiar subluminous GRB 980425, not shown in 
Fig- HI for scale compression reasons. The correlation between V 
and L is confirmed with higher significance due to the high number 
of GRBs: the linear correlation coefficient is 0.437 (P-value of 

5 X 10"^'^), the Spearman's coefficient is 0.449 (P-value of 10"^**) 
and the Kendall's coefficient is 0.302 (P-value of 2 x 10^^"). 

If we try to fit the correlation with a power law (eq. |2j, the 
result is unsatisfactory (x^/dof= 4238/549). 

logLso = m log(V7=o.45) + q (2) 

Similarly to G05 we conclude that a power-law fit is inadequate to 
describe to correlation between L and Vf =0.4,5 ■ Nevertheless, if one 
fits it with a power-law model, the best-fitting index turns out to be 
m = 0.85±0.02. This is in contrast with the value originally found 
by ROl for a sample of 13 GRBs with known redshift, m-Roi = 
3-3lJ g, and consistent with that found by G05 for a sample of 32 
GRBs with known redshift, mcos = 1.301q 44- Thus we infer that 
the correlation between L and V/=o.45 appears to be shallower than 

^ ftp://cossc.gsfc.nasa.gov/compton/data/batse/ascii_data/64ms 

^ http://cossc.gsfc.nasa.gov/batse/batseburst/sixtyfour_ms/bat_files.revamp_join 



551BNB04 0.85 ±0.02 2.45 ± 0.03 4238/549 0.70 
31 G05 l-30to;44 3.36lo 43 1167/30 0.57 



that found by ROl and this confirms what G05 found for the GRBs 
with known redshift. 

Apparently these 55 1 BNB04 GRBs seem to locate differently 
from the 31 GRBs with measured redshift (Fig.|3j. This is mildly 
suggested by the result of a K-S test, which gives D = 0.328 
with a probability of 9 x 10^^ that the two groups belong to the 
same class. If we study the scatter distribution of the 551 BNB04 
GRBs with respect to the best-fitting power-law found by G05 
for the GRBs with known redshift, with mcos = 1-301q44 
9G05 ~ 3.361q |3, it comes out that the mean residual is -0.443 (in 
Fig.|3|this corresponds to the residuals of the black circular points 
with respect to the black solid line). Hence these GRBs are on av- 
erage lO"" *'''^ ~ 0.4 times as luminous as the sample of 31 GRBs 
with measured redshift. The scatter of the 551 BNB04 GRBs is 
a = 0.7 to be compared with that of the 31 GRBs with known red- 
shift, which is 0.6. The two distributions are shown in Fig.|4] the 
black solid line shows the scatter distribution for the 551 BNB04 
GRBs, while the red dashed line shows the case of the 31 GRBs 
with measured redshift. 

The results of the best-fitting power law and the scatter for the 
two groups of GRBs (BNB04 and G05) are reported in TableQ 



5 DISCUSSION 

In agreement with the results by G05, the correlation between vari- 
ability and burst-frame 100-1000 keV peak luminosity originally 
discovered by ROl is confirmed also using a sample of 55 1 BNB04 
GRBs with a significant and consistent measure of variability, as 
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Figure 2. Vf—on^ vs. Peak Luminosity for 551 BNB04 GRBs with significant variability. Black solid lines mai'k the best-fitting power law found by ROl 
(central line) and ±1ct widths; blue dashed lines show the result that we found for 55 1 BNB04 GRBs. Red square points are GRBs with spectroscopic redshift 
taken from G05. For scale compression reasons the subluminous GRB 980425 is not shown. 
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Figure 3. vs. Peak Luminosity for 551 BNB04 GRBs with significant variability. This Figure shows the same data as Fig.|5] aside from the error 

bars, here omitted for the sake of clarity. Black solid lines mark the best-fitting power law found by G05 (central line) and ilcr widths for a sample of 32 
GRBs with known redshift (red square points). For scale compression reasons the subluminous GRB 980425 is not shown. Dashed lines show the result that 
we found for 551 BNB04 GRBs. 
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defined by ROl. Like in the case of the GRBs with known red- 
shift considered by G05, we find that the best-fitting power-law 
slope (index of m = 0.85 ± 0.02) is shallower than that found 
by ROl (mRoi = S.SlJ g). However, like G05 we find that the 
power-law description of this correlation is inadequate (x^/dof of 
4238/549). The result obtained for 551 GRBs with pseudo-redshift 
derived from the lag/peak luminosity relationship is in agreement 
with that found by G05 for a sample of 32 GRBs with spectro- 
scopic redshift and is incompatible with the best-fitting power law 
obtained by ROl. 

The comparison between the sample of 551 BNB04 GRBs and 
that of the GRBs with measured redshift considered by G05 reveals 
that the former group is on average 0.4 times as luminous as the lat- 
ter. Interestingly, despite the fact that a power law does not provide 
a good description of the correlation, the two groups show compat- 
ible power-law indices: 0.85 ± 0.02 (BNB04) vs. l.SOl 0:44 (G05). 

We conclude that the disagreement found by G05 on a sample 
of 32 GRBs with known redshift with respect to the original results 
by ROl derived from a sample of 13 GRBs with known redshift is 
confirmed by the sample of 551 GRBs with pseudo-redshifts esti- 
mated from the spectral lag/peak luminosity anticorrelation. 

A thorough discussion and possible explanations of the dis- 
crepancy between the results obtained by ROl and those by G05, 
consistent with those presented in this paper, is reported in G05. 
However, we remark that the results derived by ROl were based on 
a sample of only 13 GRBs with measured redshift available at the 
time, some of which have just limits on the peak luminosity. Re- 
markably, despite the fact that the sample of 32 GRBs with known 
redshift studied by G05 collects light curves from different space- 
craft and the sample here considered of 551 BNB04 GRBs consists 
of BATSE data only, the results are consistent with each other. 

The fact that the BNB04 GRBs are on average less luminous 
than the 31 GRBs with measured redshift could reflect the follow- 
ing possible observational bias: the more luminous GRBs are more 
likely to have measurable redshift at lower wavelengths. 

In summary, the two samples of GRBs show the same proper- 
ties in the variability/peak luminosity space: in either case the cor- 
relation is confirmed, although it turns out to be inconsistent with 
that found by ROl. The only difference, aside from the lower aver- 
age luminosity of the BNB04 sample already discussed, is a little 
higher scatter of the BNB04 GRBs around the best-fitting power- 
law: a = 0.7 to be compared with a — 0.6 of the 31 G05 GRBs. 

It is also worth mentioning that, when we selected the 551 
BNB04 GRBs out of a sample of 827 GRBs by requiring significant 
and consistent measures of variability across different energy chan- 
nels for each single GRB, a considerable fraction of the 827 GRBs 
do not match these criteria. Thus, unlike the findings by ROl, we 
find that there are GRBs whose variability as defined by ROl does 
depend on the energy channel. 

These results are n ot necessarily in contradiction with those 
bv lSchaefer et alJ <200ll) derived from a sample of BATSE GRBs. 
First of all, because the definitions of variability and peak lumi- 
nosity adopted by those authors are taken from FRROO. Regarding 
the definition of variability, from fig. 3 of FRROO we can estimate a 
power-law relation between the two measures of variability accord- 
ing to the following: Vfrr oc Vrqi , with 6 ^ 0.8 — 0.9. The best- 
fitting power law found by 'Schaefer et alj i200ll) between I/prr 
(peak luminosity as defined by FRROO) and Vprr has an index of 
msoi = 2.5 ± 1.0. Then it is: Lfrr oc V^i';^'"''^ ~ V^o\^° '' . It 
must be pointed out that the relation between the two definitions 
of variability should be investigated in detail, since the above ap- 
proximation is based on fig. 3 of FRROO, based on 8 GRBs only. 



Nonetheless, this shows that the power-law index between variabil- 
ity and peak luminosity does depend on the definitions adopted. 
Moreover, the definition of variability given by FRROO appears to 
be disputable in some points: e.g., the action of rebinning by a non- 
integer factor the light curves to report them to a fixed reference 
frame, dramatically affects the nature of Poisson counting statistics 
of the time series. Consequently, the variance that must be sub- 
tracted in the variability expression (eq. 2 in FRROO) is no more 
Poissonian and should be corrected accordingly. In general, the op- 
eration of splitting the counts integrated over a single time bin and 
use them to evaluate the variability is potentially risky, especially 
when the smoothing time-scale, proportional to the GRB duration, 
is relatively short. 



6 CONCLUSIONS 

We derived a selected sample of 551 BATSE GRBs from the 
BNB04 catalogue of GRBs assuming pseudo-redshifts based on 
the peak luminosity/spectral lag anticorrelation 1 Norris et al. 200^). 
The GRBs of this sample have been selected out of a sample of 
827 BNB04 GRBs by requiring a significant and consistent mea- 
sure of variability across the different energy channels, as defined 
by ROl. Unlike ROl, we find that not all of the 827 BNB04 GRBs 
we initially considered show a consistent measure of variability for 
different energy bands. 

We confirm the correlation between variability and peak lu- 
minosity for the subsample of 551 BNB04 GRBs. In agreement 
with the results by G05 on a sample of 32 GRBs with measured 
redshift and in contrast with the original results by ROl on a sam- 
ple of 13 GRBs with known redshift, we find that a power-law de- 
scription of the correlation is inadequate. Nonetheless, if we try 
to fit it, we obtain a power law which is significantly shallower 
(m = 0.85 ± 0.02) than that found by ROl (m = S-Sto.l) and 
consistent with that found by G05 (l.SOtoil)- 

Finally, we note that the sample of 551 BNB04 GRBs are on 
average less luminous by a factor of ~ 0.4 with respect to the GRBs 
with measured redshift. We ascribe this difference to the fact that 
the sample of GRBs with measured redshift has been biased so 
far in favour of the more gamma-ray luminous GRBs, which on 
average had high fluences and then for which, consequently, more 
precise localisations were possible. This is probably true for the 
very first GRBs with measured redshift. We expect that Swift will 
clarify this issue as soon as it will discover a number of GRBs with 
low peak luminosity and measurable redshift. 
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